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Abstract 

In this paper we prove that Bose gas with weak pair interaction is 
non ergodic system. In order to prove this fact we consider the diver- 
gences in some nonequihbrium diagram technique. These divergences 
are analogous to the divergences in the kinetic equations discovered by 
Cohen and Dorfman. We develop the general theory of renormaliza- 
tion of such divergences and illustrate it with some simple examples. 
The fact that the system is non ergodic leads to the following conse- 
quence: to prove that the system tends to the thermal equilibrium we 
should take into account its behavior on its boundary. In this paper 
we illustrate this thesis with the Bogoliubov derivation of the kinetic 
equations. 
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1 Introduction 

In this paper we study divergences in some nonequilibrium diagram technique 
which is analogous to the Keldysh diagram technique. It is more or less 
evident that these divergences are the same as the divergences in the kinetic 
equations discovered by Cohen and Dorfman [1,2]. We develop the general 
theory of rcnormalization of such divergences analogously to the Bogoliubov 
— Parasiuk /^-operation method [3-5]. Our main result can be formulated 
as follows: for a wide class of Bose systems in the sense of formal power 
series on coupling constant there exists non-Gibbs functional (•) such that 
the correlators 

(^±(t,xi)...^±(t,a;„)) (1) 

are translation invariant, do not depend on t and satisfy the weak cluster 
property. Here are the fields operators and the weak cluster property 
means the following 

lim / {^^{t, xi + 5ieia)...'^^{t, Xn + 5neia))f{xi, Xn)d^xi...(fxn 

\a\^oo J 

= j {^^{t,Xi,)...^^{t,Xi,)){^^{t,Xi,)...^^{t,X,^)) 

Xf{xi,...,Xn)d^Xi...d^Xn, (2) 

there 6i G {1,0}, i = l,2...n and 

ii <i2 < ■■■ < ik, 

ik+l < ik+2 < ■■■ < in, 

{^l, ^2, ik} = {i = 1, 2...n\5i = 0} ^ 0, 
{ik+i, ik+2: in} = = 2...n\5i = 1} ^ 0. (3) 

/(xi, ...,Xn) is a test function, ei is a unit vector parallel to the a;-axis. This 
statement is a simple consequence of the theorem from the section 6. 

Let us prove that the existence of such functionals implies non-ergodic 
property of the system. We consider the problem only on classical level. 
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The accurate consideration for the quantum case can be found in section 
10. Suppose that our system is ergodic, i.e. there are no first integrals of 
the system except energy. Then, the distribution function depends only on 
energy. We can represent the distribution function f{E) as follows: 

= J] c„5(£; (4) 

where the sum can be continuous (integral). Let 1 be some enough large 
but finite subsystem of our system. Let 2 be a subsystem obtained from 

— * 

1 by translation on the vector I of sufficiently large length parallel to the 
X-axis. Let 12 be a union of the subsystems 1 and 2. Let pi, p2 and pi2 
be distribution functions for the subsystems 1, 2 and 12 respectively. Let 
Fi, r2 and be points of the phase spaces for the subsystem 1, 2 and 12 
respectively. By the same method as the method used for the derivation of 
the Gibbs distribution we find: 

Cada^ ^ , C?a > OVa (5) 

in the obvious notation. But the weak cluster property implies that 

Pl2 = PlP2- (6) 

Therefore all the coefficients Cq are equal to zero except one. We find that 

f{E) = c6{E - Eo) (7) 

for some constants c and Eq. So each finite subsystem of our system can be 
described by Gibbs formula and we obtain a contradiction. 

Non-ergodic property means that there is no thermalization in infinite 
Bose-gas system. 

This fact implies to prove that the system tends to thermal equilibrium 
we should take into account the behavior of the system on its boundary. 
Indeed if a system has no boundary the system is infinite. 

To illustrate this fact we will study Bogoliubov derivation of kinetic equa- 
tions [6] . When one derives BBGKI-chain one neglects some boundary terms. 
If one takes into account this boundary terms and uses the Bogoliubov 
method of derivation of the kinetic equations one finds that these bound- 
ary terms compensate the scattering integral. 
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I think that the dependence of behavior of the system of boundary can 
be observed for small systems such as nanosystems or biological systems. 

Note that our main result is closely related with so-called the Prigogin 
hypothesis which states that the infinite dimensional Liuville dynamics can 
not be derived from the Hamilton dynamics. The Prigogin hypothesis is 
proven in [7] 

The paper is organized as follows. In section 2 we introduce the notion of 
the algebra of canonical commutative relations and develop an useful repre- 
sentation for some class of the states on this algebra. In section 3 we describe 
the von Neumann dynamics for the states. In section 4 we describe an useful 
representation for the von Neumann dynamics — the dynamics of correla- 
tions. In section 5 we describe the decomposition of the kinetic evolution 
operator by so-called trees of correlatios. In section 6 we describe the gen- 
eral form of the counterterms which subtract the divergences in the nonequi- 
librium perturbation theory. In section 7 we describe so-called Priedrichs 
diagrams. In section 8 we describe the Bogoliubov — Parasiuk prescriptions 
and formulate our main theorem. In section 9 we prove our main theorem. In 
section 10 we derive the non-ergodic property of Bose gas with weak pair in- 
teraction from our main result. In section 11 we consider one example related 
to our general theory. In section 12 we reconsider the Bogoliubov derivation 
of the Boltzmann equation. This example illustrates the main thesis of this 
paper: to prove that the system tends to the thermal equilibrium one has to 
take into account its behavior on its boundary. Section 13 is a conclusion. 

2 The Algebra of Canonical Commutative Re- 
lations 

Let S'(R'^) be a Schwatrz space of test functions (infinitely-differentiable func- 
tions decaying at infinity faster than any inverse polynomial with all its 
derivatives). The algebra of canonical commutative relations C is an uni- 
tal algebra generated by symbols a+(/) and a(/) / e satisfying the 

following canonical commutative relations: 

a) a^{f) is a linear functional of /, 

b) a{f) is an antilinear functional of /, 

[a(/),a(^)] = [a+(/),a+(/)] = 0, 
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Hf),a^g)]^{f,g), (8) 
where {f,g) is a standard scalar product in L^(M^), 

{f,g):^ j d'xr(x)g(x). (9) 

Let po be a Gauss state on C defined by the foUowing correlator 

Poia+{k)a+{k')) = po{a{k)a{k')) = 0, 

PQ{a+{k)a{k')) = n{k)S{k - k'), (10) 

where n{k) is a real-valued function from the Schwartz space. In the case 
then 

g-/3(a;(fc)-M) 

^(^) = 1 _ e-/3(^(fe)-.) ' (11) 

where // e IR, // < 0, po is called the Plank state. Here uj{k) — 

Let C be a space of linear functionals on C, and C'|_ be a set of all states 
on C. Let us make the GNS construction corresponding to the algebra C and 
the Gauss state po- We obtain the set {7i,D,^,)) consisting of the Hilbert 
space Tt, the dense linear subspace D in H, the representation "of C by means 
of the linear operators from D to D, and the cychc vector ) e D, i.e. the 
vector such that C) — D. This set satisfies the following condition: Wa E C 
(a) = po(a)- Below we will omit the symbol", i.e. we will write a instead of 
a. 

Let us introduce the field operators: 

*(^) = / ^"'Mk)dk, 

[271)2 J 

^+(a;) = [ e-'^''a+{k)dk. (12) 

(27r)2 J 

We say that the state p on C satisfies the weak cluster property if 
lim / {^^{t, xi + 5ieia)...*=^(i, x„ + (5„eia))/(xi, Xn)d^Xi...(fxn 

a— »oo J 

= j (vI/±(t,a:,J...*±(t,x,J)(*±(i,x,J...*±(i,x,J) 

Xf{xi,...,Xn)d^Xi...d^Xn, (13) 
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where Si G {1,0}, i = l,2...n and 

ii <i2 < ... < ik, 

ik+l < ik+2 < ... < in, 

{ii,i2, ■■■,ik} = {i = l,2...n\6i = 0} 7^ 0, 
{ik+i,ik+2, in} = {i = l, 2...n\Si = 1} 7^ 0. (14) 

f{xi, Xn) is a test function, ei is an unit vector parallel to the a;-axis. 
Definition The vector of the form 

j v{pi, ...,Pn)a^{pi)...a^{Pn))d^Pi---d^Pn, 

v{pu...,Pn)eS(R^^). (15) 

is called a finite vector. The finite linear combination of the vectors of the 
form f|T5l) is also called a finite vector. 

Let f{xi, Xklvi, yi\vi, Vmlwi, Wn) be a function of the form 

/(Xi, ...,Xk\yi, --^yilVl, ...,Vm\Wi, ...,Wn) 

= g{xi, ...,Xk\yi, ■■■,yi\vi, ...,Vm\wi, w„) 

k m I n 

x5(^Xi-^Wj-^l// + ^t;3), (16) 

i=l j=l f=l g=l 

where g is a function from Schwartz space. 
Consider the following functional on C 

/k I m n 

Y[ dxi JJ dxj Wdvf^ dwg 

i=l i=l /=1 9=1 
Xf{xi, ...,Xn\yi, ■■■,yi\Vl, ...,Vm\Wi, ...,Wn) 
, s 

xpo(: a[xi)...a[xn) {y\) . . .a^ {yi) : A : a{vi)...a{vrn) a^{wi)...a^{wn) :). (17) 
Here the symbol 

: (.T~^~T-) : (18) 

means that when one transforms the previous expression to the normal 
form according to the Gauss property of po one must neglect all correla- 
tors po('3'^(xi)a^(x„) such that a^(xi) and a^(x„) both do not come from 
A. 
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Let C be a subspace in C spanned on the functionals just defined. 
Now let us introduce an useful method for the representation of the states 
just defined. 

Let C2 = C+(8>C_, where C+ and C_ are the algebras of canonical commuta- 
tive relations. The algebras C± are generated by the generators a±{k), a^(fc) 
respectively satisfying the following relations: 

[aXik),aXik')] = [a4k),a4k')] = 0, 
[atik),atik')] = [a^{k), a^{k')] = 0, 

[a+{k),aX{k')] = 6{k-k'), 

[a-{k),at{k')] = 6{k-k'), 

[aX{k),a^{k)] = 0. (19) 

Here we put by definition := a±. Let us consider the following Gauss 
functional p'^ on C2 defined by its two-point correlator 

p'o{a^ik)a^ik')) = po{a^ik)a^ik')), 
p'oiaXik)aXik')) = poia^ik')a^ik)), 
p',{aX{k)aZ{k')) = p',{a-{k)at{k')) = 0, 
p'Q{al{k)aZ{k')) = n{k)6{k - k'), 
Po{aX{k)at{k')) = (1 + n{k))6{k - k'). (20) 

One can prove that the functional p'q is a state. 

Let us make the GNS construction corresponding to the state Po and the 
algebra €2- We obtain the set {Ti' , D,'^, )) consisting of the Hilbert space H', 
the dense linear subspace D in Ti.', the representation " of C2 by means of the 
linear operators from D to D, and the cyclic vector ) G D, i.e. the vector such 
that C) = D. This set satisfies the following condition: Va G C2 (a) = p'q^o). 
Below we will omit the symbol i.e. we will write a instead of a. 

Now we can rewrite the functional, defined in f|T71) pf as follows 

Pf{A) = (A'Sf), (21) 

where A' is an element of C2 such that it contains only the operators a_ , al 
and can be represented through a-,a1i in the same way as A can be repre- 
sented through a, a"*". Sf is an element of C2 of the form 

/k I m n 

Y\ W dxj Y\ dvf JJ^ dwg 
*=i i=i /=i 9=1 
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X/(X1, ...,Xn\yi, --^yiW, ...,Vm\Wi, ...,Wn) 

X : aX{xi)...aXiXn)a+{yi)...a+{yi)a^{vi)...a_{v^)at{wi)...at{wn) : • (22) 

Here the symbol : ... : is a normal ordering with respect to the state pg- 
Denote by D' the space dual to D. We just construct the injection from 

C into D'. Denote its image by 7i'. 

By definition the space C" is a space of all functionals on C which can be 

represented as finite linear combinations of the following functionals 

p{A) = {A':Sf,...Sf,^:). (23) 

Here A' is an element of C2 such that it contains only the operators a_, and 
can be represented through a_, al in the same way as A can be represented 
through a, a"*" and Sf. are the elements of the form fl22l) . Denote by Ti." the 
subspace in D' spanned on the vectors : Sfj^...Sf„ :) (in obvious sense). 
There exists an involution -k on 7i' defined by the following formula: 

/k I m n 

Y\ '^^i Y\. '^^j n '^^f n '^^3 
i=l j=l f=l 9=1 

X/(xi, ...,Xn\yi, ---^yilVl, ...,Vm\Wi, W„) 

X : aX{xi)...aX{xn)a+{yi)...a+{yi)a^{vi)...a^{vm)at{wi)...at{wn) :)}* 

k I m n 

W dxi JJ^ dxj JJ^ dvf JJ^ dwg 

i=l i=l /=! 9=1 

x/*(a;i, ...,Xn\yi, ■■■,yi\vi, ...,v„,\wi, ...,Wn) 
X : atixi)...at{Xn)a.{yi)...a_{yi)a+{vi)...a+{v.ra)aX{wi)...aXiWn) :). (24) 

We define the involution -k on Hom(if ', H') by the following equation: 

{a\f)y = a*m\ (25) 

where a G Hom(^', if') and |/) G H'. 

We define also the involution -k on by the following equation: 

/k I m n 

Y\ dXi JJ^ dxj JJ^ (ii;/ JJ^ rfty^ 
«=i i=i /=i 9=1 
x/(xi, ...,w„) 
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X : aX{xi)...aX{xn)a+{yi)...a+{yi)a^{vi)...a^{vm)at{wi)...at{wn) :}* 

k I m n 

Y\ Y\. '^■^■^ n ■i^ n 

«=i i=i /=i 9=1 
x/*(a;i, ...,a;„||/i, ...^yilvi, ...,v^\wi, 

X : atixi)...at{xn)a_{yi)...a_{yi)a+{vi)...a+{vra)aX{wi)...aX{wn) :, (26) 

where /(xi, ...,Xk\yi, t/ilfi, ...,^^1^1, w„) is a test function of its argu- 
ments. Note that the involution on Hom(i7', H') extends the involution on 
C^. We say that the element a e C2 is real if a* = a. The involution on H" 
can be defined by a similar way. 



3 The von Neumann Dynamics 

Suppose that our system is described by the following Hamiltonian 

H^Ho + XV, (27) 

where 

ifo = y d?k{ixi{k) - ii)a^{k)a{k) and 

xd{pi +p2-qi- q2)a^{pi)a^{p2)a{qi)a{q2). (28) 

Here the kernel v{j)i,p2\qi, Q'2) belongs to the Schwartz space of test functions. 
To point out the fact that H is represented through the operators a"*", a~ we 
will write H{a'^, a~). 

The von Neumann dynamics takes place in the space H" and is defined 
by the following differential equation: 

(29) 

where the von Neumann operator has the form 

-iH{at,aZ)+iH\aX,a~), (30) 
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where we put by definition: 

/n m 
Y\dpiY\_dqjv{pi, ...,pn\qu ...,qm) ■ a+(pi)...a+(pn)a(gi)...a(?„) :)^ 
i=i j=i 

/ n m 
YldpiY[dqjv(pi, ...,pn\qi, ...,qmT ■ a'^ipi)--a'^iPn)a{qi)--a{Qn) ■ -(31) 
i=i j=i 

Let us divide the von Neumann operator into the free operator C and the 
interaction dnt, C — jCo + ^J^int, where 

£,0 = -iHo{a1l, aZ) + ii/o(a|, a+), 
Cint = -iHint{at, aZ) + iHj^^{aX, a~). (32) 

Note that the operators Cq and Ci are real (with respect the involution 

Let us introduce kinetic evolution operator (in the interaction represen- 
tation) 

C/(t",i') = e-^°*"e^(*"-*')e^°*'. (33) 

After differentiating with respect to t we find the differential equation for 
U{t,t'). 

jUit,t')^jCintmit,t'), (34) 

where 

jCintit) = e-^°*A„te^°*. (35) 

So the state )p under consideration in the space Ti" in the interaction repre- 
sentation has the form 

o 

)p = rexp( J Cint{t)dt)), (36) 

— oo 

where T is the time-ordering operator. 

Note that we have a hnear map from H" into C . It is easy to see that 
the von Neumann dynamics is in agreement with the Heizenberg dynamics 
in C. 
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4 Dynamics of Correlations 



Let us construct some new representation of the von Neumann dynamics 
useful for the renormahzation program. This representation is called the 
dynamics of correlations. The ideas of the dynamics of correlations belongs 
to I. Prigogin [8]. The dynamics of correlations takes place in the space 

oo 

K:=0sym(g)"7i:'. (37) 



Now let us describe how the operators Cq and C%^^ act in the space Tic- 

Let us define the actions of operators £q and C'l^^ which are corresponds 
to the operators Cq and Cint- 

By definition all the spaces (8)"7Y' are invariant under the actions of op- 
erators Cq. Note that the space li! is invariant under the action of operator 
Cq. Let us denote the restriction of Cq to the space Hf by the symbol C'q. 
By definition the restriction of Cq to the each subspace sym (g)" Ti' of Tic has 
the form 

£() 1 (8) ... (g) 1 + 1 £o (8) ... 1 + ... + 1 (8) 1 ... (8) £(). (38) 

Now let us define C'l^^. Let |/) e Tic, belongs to the subspace (8)"7Y' and has 
the form: 

m 

(39) 

i=0 

where /j has the form 

/k I m n 

n n n n 
i=l j=l f=l g=l 

Xf{xi, Xn\yi, yi\vi, Vm\wi, Wn) 

X : aX{xi)...aX{xn)a+{yi)...a+{yi)a-{vi)...a-{vrn)a1:iwi)...at{wn) : . (40) 
By definition, 

^tt\f)-0 (41) 
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if / > n. Let us consider the following vector in TY" 



m n 



H^Wn^- (42) 

i=i j=i 

m n 

Let us transform the expression Cint X] ■ 11 /j ■ ^^e normal form. Let us 

i=i j=i 

denote by hi the sum of all the terms in the previous expression such that 
exactly I operators /j couple with Cint- We find that hi) has the following 
form 

k 

hi) = J2-- 9\-9n-i+i ■■) (43) 

i=l 

for some k. Here gl has the form of right hand side of fj40l) . Now let us 
consider the following vector 

k 

\f)i = sym J] : gi :) ® : gl^_i^^ :), (44) 

i=i 

where we define symmetrization operator as follows 

sym(/i (g) ... (g) /„) 

= ;^E^-®-® (45) 

— the group of permutation of n elements.) Put by definition 

^itlf) = \f)i- (46) 

One can prove that this definition is correct. Analogously, in the following 
expression 



m n 



An^E^n^^) (47) 

i=l j=l 

let us keep only the terms such that Cint does not couple with any of /j. Let 
us write the sum of such terms as follows 

E^n^5o- (48) 

i=i j=i 
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Here /i*) has the form of right hand side of fl40p . Let \h) be a vector in 
sym®"^^7i' defined as follows 

/ n+i 

\h)=symY,®--h]:). (49) 

i=\ j = i 

Put by definition 

^^t\f) = \h). (50) 

We have the evident linear map F : Ti^ — > H" which assigns to each 
vector sym : /i :) (g) ...(g) : :) the vector : /i.../n :)• Denote by If^ the evo- 
lution operator in interaction representation in the dynamics of correlation. 
The following statement describes the relation between the von Neumann 
dynamics and the dynamics of correlations. 

Statement. The following relation holds: 

FoU'{t',t") = U{t',t")oF. (51) 

5 The tree of correlations 

The useful representation of dynamics in Tic is a decomposition by so called 
trees of correlations. 

Definition. A graph is a triple T = {V, R, /), where V, R are finite sets 
called the set of vertices and set of lines respectively and / is a map: 

h: R^V^^^ UV X {+}UV X {-}, (52) 

where V^'^^ is a set of all disordered pairs {vi, V2), Vi, V2 E V such that Vi 7^ V2- 
If (fi,f2) = f{r) for some r G -R we say that the vertices Vi and V2 are 
connected by a line r. If /(r) = (fi,t>2), vi,V2 G V we say that the line r is 
internal. 

Remark. We use this unusual definition of graphs only in purpose of 
this section to simplify our notations. 

Definition. The graph F is called connected graph if for two any vertices 
V, v' there exists a sequence of vertices v = vq, vi, ...,Vn = v' such that = 
0, n — 1 the vertices Vi and fj+i are connected by some line. 

By definition we say that the line r is an internal line if /(r) = (^1,^2) 
for some vertices Vi and V2- 
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For each graph T we define its connected components by the obvious way. 

Definition. We say that the graph F is a tree or an acychc graph if the 
number of its connected components increases after removing an arbitrary 
hne. 

Definition. The elements of the set f~^{V x { — }) we call the shoots. 
Put by definition R^h = f~^{V x {-}). The elements of the set f-'^{V x {+}) 
we call the roots. Put by definition Rroot = /"^{V x {+})• 

Definition. Directed tree is a triple (T, $„, ^sh), where T is a tree and 

and are the following maps: 

^sh-Rsh^{l,'^,-ARsh}- (53) 

Definition. We will consider the following two directed trees (T, ^gh) 
and (T", $'^^) as identical if we can identify the sets of lines R and R' of 

T and T' respectively and identify the sets of vertices V and V of T and T' 
respectively such that after these identification the trees T and T' become 
the same, the functions and $^ become the same and the functions 
and ^'gf^ become the same. 

Denote by r(T) the number of roots of T and by s{T) the number of 
shoots of T. Below, we will denote each directed tree (T, by the 

same symbol T as a tree omitting the reference to and write simply 

tree instead of the directed tree. 

We say that the connected directed tree T is right if there exists exactly 
one line from f~^{V x {+}). 

We say that the tree T is right if each its connected component is right. 

The vertex v of the tree T is called a root vertex if (v, +) G f^^{R). 

To point out the fact that some object A corresponds to a tree T we 
win often write At- For example we will write T = (Vr, Rt, It) instead of 
T^{V,RJ). 

Definition. For each connected right tree T there exists an essential 
partial ordering on the set of its vertices. Let us describe it by induction on 
the number of its vertices. Suppose that we have defined this relation for 
all right trees such that the number of their vertices is less or equal than 
n — 1. Let T be a right tree such that the number of its vertices is equal to 
n. Let Vmax be a root vertex of T. Put by definition that the vertex v^ax is a 
maximal vertex. Let be all of its children i.e. the vertices connected 

with Vmax by lines. By definition each vertex vi < Vmax,i = l,---,k. We 
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can consider the vertices Vi,...,Vk as a root vertices of some directed trees 
Tj, i = 1, k. By definition the set of vertices of Tj consists of all vertices v 
which can be connected with Vi by some path v = v'l, = fj such that 

Vmax 7^ for all J = 1,...,/. The incident relations on Tj are induced by 
incident relations on T. Put by definition that i,j — l,...,k, i ^ j 

and for any two vertices v[ G Tj and V2 G Tj v[ ^iv^- If G Tj for some 
Tj we put v'^ ^ f 2 in T if and only if v'^ ^ v'2 in the sense of ordering on Tj. 
We put also v < Vmax for every vertex v 7^ v^ax- These relations are enough 
to define the partial ordering on T. 

If the tree T has several connected components we define a partial order- 
ing at each its connected components as previously and put Vi V2 if Vi and 
V2 do not belongs to the same connected component of T. 

Below without loss of generality we suppose that for each tree of cor- 
relation T and its line r the pair (^1,^2) = /(r) satisfies to the inequality 

Vi > V2. 

Definition. The tree of correlations C is a triple C = (T, (p, r), where T 
is a directed tree, r is a map from R \ Rgh to := {x G ]R|x > 0}: 

{r{r)UR = r{r), (54) 
and </? is a map which assigns to each vertex v of T an element 

ip{v) eRom{(^n','H') (55) 

r—*v 

of a space of linear maps from 'H' to H'. 

r — >v 

In H' the tensor product is taken over all lines r such that r ^ v. 

(r— 

Let V be a vertex of the tree T. If /(r) = {v',v) for some vertex v' or 
/(r) = [v, +) we say that the line comes from the vertex v and write r -s— v. 
If /(r) = (■;;, v') for some vertex v' or /(r) = {v, — ) we say that the line comes 
into the vertex v and write r — > v. 

Definition. Let (T, 99, r) be a tree of correlations such that for each 
vertex v (f{v) — CfJ^^, where ly is a number of lines coming into v. We call 
this tree the von Neumann tree and denote it by Tf. We also say that (f is 
a von Neumann vertex function. 
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Definition. To each tree of correlations {T, cp, r) we assign an element 



t/*,^(f)eHom(07Y',07i:') 



(56) 



by the following way: 

If T is disconnected then 



®{uhT,c^m (8) fi}- 



(57) 



CT ieR,h{CT) 



Here the number of connected components of T is equal to n, and connected 
components of T are denoted by CT. Ccp and Cr are the restrictions of (p 
and r to the sets of vertices and lines of CT respectively. Rsh{CT) is a set 
of shoots of CT. Now let T be a connected tree. To define 



by induction it is enough to consider the following two cases, 
case 1). The tree T has no shoots. 

a) Suppose that the tree T has more than one vertex. Let Vmin be some 
minimal vertex of T and vq be a vertex such that an unique line tq comes 
from Vjnin into Vq. Let T' be a tree obtained from T by removing the vertex 
fmin of T. Let -f be a restriction of r to \ {tq}. Let ip' be a function, 
defined on y \ {vmin} as follows: ip'{v) — (p{v) ii v ^ vq and 



UtJ^) (g) fr 



(58) 



reRsh 



f'iVo) (g) fr^ ^{Vo) (g) h, 



(59) 



where 



hr = fr a r To, and 



' ''To — ^ r\ '^mm ) • 



(60) 



Put by definition 



(61) 
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b) The tree T has only one vertex Vmin- Then 

UT,^{r) = e-(*-^)^V(^^m^n). (62) 

Case 2.) The tree T has a shoot tq coming into the vertex vq. In this case 
instead of the tree {T,(f,T) we consider the tree {T',lp',t'), where the tree 
T' has the same vertices as T, the set of hnes of T is obtained by removing 
the hne Tq from the set of hnes of T', the function r' is a restriction of the 
function r to the set of hnes of T' and the function </?' is defined as follows: 

(p'{v) = (p{v), iff 7^ fo and 
v'i'^o) hr = (p{vo) g-r, where 

1 — >?;o; r^ro i — >vo; 

Qr — hr, ii r ^ To, and 

5r = e^°(*-*^)^. (63) 

Here we put 

tr^J2^r', (64) 

where the sum is taken over all lines r' which forms decreasing way coming 
from + to Vq. Put by definition 

U'Tjr)\f)---U'^',A^')\f), (65) 

where 

1/')= (g) (66) 

re(i?sh)T' 

Let (T, (/?, r) be some tree of correlations. We can identify the tensor 
product 



n'r (67) 

reRsh 



with 

sh{T) 



n' (68) 

i=l 



17 



and the tensor product 



(g) K (69) 

reRroot 

with 

r(T) 

^n.' (70) 

Using these identifications let us consider an operator V^^(t) : 'W 'W 
defined by the foUowing formula 

T4_^ = symoC/*,^^oP,^(r), (71) 

sh{T) ^ 

where Psh{T) is a projection of Tic to sym H.' . 

1=1 

Remark. If (T, 9?, r) is a von Neumann tree of correlations then we will 
shortly denote the operators f^(V<^) ^^^d y^Ttp) ^^<i ^ respectively. 

The following theorem holds: 

Theorem. The following representation for the evolution operators holds 
(in the sense of formal power series on couphng constant A). 

U'{t',t") = J2— / V^{f)dr. (72) 

WeRsh t-tr>t" 

Here ut is a number of vertices of the directed tree T. 



6 The general theory of renormahzation of 

U{t,-oo)) 

In the present section we by using the decomposition of correlations dynamics 
by trees describe the general structure of counterterms of U{t,— 00)), which 
subtract the divergences from U{t, —00)). We will prove in the section 10 
below that there exist divergences in the theory. Note that the structure of 
i?-operation for the processes at large times for some class of systems has 
been considered at [9] 
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Let T be a tree. Let us give a definition of its right subtree. 

Definition. Let fi, f„ be vertices of T such that Vi, j = 1, ...,n, i ^ j 
Vi ^ Vj. Let us define subtree T„^,...,i,„. By definition the set of vertices 
^Ty^,...,v„ of consists of all vertices v such that v < Vi for some i = 

l,...,n. 

The set Rt of all lines of the tree T„. ^ consists of all lines r of Rt 
such that h{r) = {v",v') and v'.v" ^ Vi for some i = l,...,n. The incident 
relations on Ty^^,,,^vn ^.re induced by the incident relations of T except the 
following point: if the line r comes from the vertex v into Vi, i = 1, ■■■,n we 
put fT^^^ v„}i^) ~ ('^' +)• ^^i^ ^i^® r is a root and the vertex v is 

a root vertex of the tree Vrr , . The tree TL, „ is called a right subtree 

^ {vi,...,Vn} 01,...,Un O 

of T. 

The Bogoliubov — Parasiuk renormalization prescription. Let 

us define the following operator: 

Wr,{t)^ (g) Zr,r,{t), (73) 
reRroot{T) 

where by definition, 

Zr,ro{t) = 1, if 7^ J^o, and 

Zr,r,{t) ^ e-'-^K (74) 

We say that the amplitudes {At^^p) are time — translation invariant ampli- 
tudes if for each tree T and for each its root fine Tq 

WrM^T,^ = ^T,^- (75) 

For each set of amplitudes >lr,(^ put by definition: 

At,^) ^ F o At,^, (76) 

where T is an arbitrary tree without shoots. 
Now let us formulate our main result. 

Theorem. There exists a procedure called renormalization which to each 
tree T without shoots assign the amplitudes At,<^ satisfying to the following 
properties a)-f): 

a) If the tree T is not connected and {CT} is a set of its connected 
components, while {dp} is a set of its restriction of </? to CT 

At,^ = 0Act,c^ (77) 



19 



in obvious notations. 

b) The amplitudes A^ ,^ are real i.e. 

(At,^)* = At,^. (78) 

c) The amplitude A^ ,^ satisfies the property of time-translation invari- 
ance. 

It has been proven that 

U{t,-oo)) = J2— drU'^r)). (79) 

^ J 

In the last formula the summation is taken over all trees T without shoots. 
Let T be a tree without shoots and T' be a right subtree of T in the 
described before sense. Let us define the amplitude 

^t',^^U't,^{t). (80) 

Let by definition T\T' be a tree obtained by removing from the set Vt all 
the vertices of T' and from the set Rt all the internal lines of T'. In (1791) r 
is a map from Rt\t' into M^. 

We can consider the amplitude U^\rp, as a map 



H'^ H'. (81) 

(■^T\T' )s/i C-^T\T' )root 



By using this identification we simply put 



= U'T\T'Jr)AT',^. (82) 

Now let us define the renormalized amplitudes, by means of the counterterms 
At by the following formula: 

(i?A?7)(t,-oo) 

\n.T f 

= E ^TT 5Z / ^ UTir)df. (83) 
d) The renormalized amplitudes {RtJJ){t, — oo)) are finite. 
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f) Let T be an arbitrary connected tree without shoots. Consider the 
following element of H': 

a-=Y.Yl f^^'^ U^ir))df. (84) 
We can represent the element a as follows: 

oo „ 

a= / Wm{xi, ...,Xk\yi, ...,yi\vi, ...,Vf\wi, ...,Wg) 

k,l,f,g=0'^ 

km Im fm Qm 

: ^^aX{xi)dxi'^a+{yi)dyiY\^a-{vi)dvi'^at{wi)dwi :). (85) 

i=l i=l i=l i=l 

Let Wk,ij^g{zi, .... Zn) (n = km + Im + fm + Qm) bc a Fourier transform of 
Wk,i,f,g{xi, ...,Xk\yi, ■■■,yi\vi, ...,Vf\wi, ...,Wg). Then 

J dzi, dZnWk,ij,g{zi + s(l)eia, + s(n)eia)/(^i, ^„), (86) 

tends to zero as a as a ^ +oo. Here s{i) are the numbers from {0, 1} and 
there exist numbers j, i^j = 1, ...,n such that s{i) = 0, s{j) = 1 for some 
i,j = l,...,n. f{zi,...,Zn) is a test function, ci is a unit vector parallel to 
the X-axis. 

Remark. The property d) imphes the weak cluster property of the func- 
tional {RU){0,-oo). 

This theorem is a simple consequence of the theorem-construction from 
the section 8. 

The renormalizcd amplitudes satisfy to the following properties: 
Property 1. For each t eM. 

(i?At/)(t,^00)) 

= e-^"*(i?AC/)(0,-oo)). (87) 

This property simply follows from the definition of {Rj'yU){t, — oo)) and means 
that the functional {RAU){t, — oo)) is a stationary state. 
Property 2. 



{RAU){t,-oo)) = t/(i,0)(i?At/)(0,-oo)). (88) 
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This property follows from the following representation of {R\U){t, — oo)). 



{RAU){t,-oo)) = U{t,-oo)I), 



(89) 



where 



T ^ 



(90) 



and the sum in the last formula is taken over all von Neumann trees without 
shoots. Property 2 means that the functional {R\U){t, — oo)) satisfies to the 
von Neumann dynamics. 

Remark. The existence of the stationary translation invariant functional 
satisfying to the weak cluster property follows from the previous theorem and 
the properties 1, 2. 

7 The Friedrichs diagrams 

Now let us start to give a constructive description of the counterterms 
such that the amplitude R{U){t, —oo)) is finite, and the counterterms 
satisfy the properties a) — f) from the previous section. 

At first we represent fZ-f ^(r), where T is some tree without shoots, as a 
sum taken over all so-called Friedrichs graphs $ concerned with T. 

Definition. A Friedrichs graph $t concerned with the directed tree T 
without shoots is a set (V ,R,Or, f^, f~ ,g), where F is a union of the set 
of vertices of T and the set {©}• Recall that there is a partial order on Vt- 
We define a partial order on the set V if we put Vf G © > f . and 

are the maps f^,f~ : R ^ V such that f~^{r) > f~{r). Or is a map 
i? ^ {+, — } called an orientation, gi is a function which to each pair (w,r), 
V e Vt, r E R such that f'^{r) — v or f~{r) — v assigns + or — . The graph 
{V, R, Or, /"•", /~, g) must satisfy the property: if we consider © as a vertex, 
the obtained graph is connected. 

If f~^{r) = V we write r — > and if /~(r) = v we write r <— 

If we want to point out that the object B concerned with the graph $ we 
will write B^. For example we will write V$ and for the sets of vertices 
and lines of $ respectively. 

At the picture we will represent the elements of V by points and the 
element © by ©. We will represent the elements of R by lines. The line 
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r connects the vertices f'^{r) and /^(r) at the picture. We will represent 
orientation Or{r) by arrow on r. If Or{r) = + the arrow is directed from 
f~{r) to f~^{r). If Or(r) = — the arrow is directed from f~^{r) to f~{r). To 
represent the map g : {r,v) —>■ {+, — } we will draw the symbol g{{r,v)) (+ 
or — ) near each shoot (r, v). At the picture a shoot (r, v) is a small segment 
of the line r near v. 

Definition. The Friedrichs diagram F is a set (T, $, (/?, /i), where T is a 
tree, $ is a Friedrichs graph, is a map which assigns to each vertex v of T 
a function of momenta {pr\r G R<s>} of the form 

ji 

Si i=l 

where is a test function of momenta coming into (from) the vertex v. 
{Si}^=i is a decomposition of the set of shoots of v into riy of disjunctive 
nonempty sets S^, p], ...,pf are momenta corresponding to the shoots from 
{Si}, h is a function which assigns to each pair v E V, r E R such that 
/+(r) > V > f-{r) a real positive number h{v,r). 

It will be clear that it is enough to consider only the diagrams F such 
that for each its vertex v and set Si e {Si}il^ there exists a line r such that 
{r,f-ir))eSi. 

To each Friedrichs diagram F = (T, (/?) we assign an element of H'^ of 
the form 

U{T,'S>,ip) (^) 

= J ...dpe,t...U'r{...Pr.^,...) (92) 
X : ...atiPr^J... :). (93) 

Here pr^^^ are momenta of external lines, i.e. such lines r that f^{r) = ©. We 
choose the lower index of a±(prext) by the following rule. Let v be a vertex 
such that f~{rext) — v. li g{{r,v)) — -\- we choose + as a lower index, and 
if g{{r,v)) = — we choose — as a lower index. We choose the upper index 
of a±(Prea:t) ^y thc following rule. If the lower index of ci%{Pr^.^t) {"} then 
the upper index is equal + if the corresponding line comes from the vertex v 
and this index is equal — if the corresponding line comes into the vertex v. 
If the lower index of a±{Prext) is {+} then the upper index is equal — if the 
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corresponding line comes from the vertex v and this index is equal + if the 
corresponding line comes into the vertex v. 

Now let us describe the amplitude t/p(...pea;t---)- By definition we have 



C/°(f)(...Pe.t...) 




n G(Or(r), giir, /+(r))), ^((r, /-(r))))(p). 



(94) 



Let us describe the elements of this formula. i?r is a set of all lines of 
diagram F. Symbol r ^ v denotes that the line r comes into (from) the 
vertex v. In the expression 



we take the upper sign + if the line r comes into the vertex v and we take lower 
sign — in the opposite case. The symbol Rt denotes the set of lines of the tree 
T from the triple (T, $, ip) and symbol tt means the line from Rt- The sym- 
bol Vr denotes the set of all vertices v such that f^{r) > v > f^{r). The sym- 
bol {RT)r denotes the set of all lines tt of Rt such that the increasing path 
coming from /~(r) into /^(r) contains tt- G{Or{r), g{f^{r)), g{f^{r))){p) 
is a factor defined as follows 



Below we will simply write Gr{p) instead of G{Or{r), g{f'^{r)), g{f~{r))){p). 

It is evident that we can represent U^{t) as a sum taken over some 
Friedrichs diagrams F corresponding to the tree T of the quantities U^is). 

Now let us define the quotient diagrams. 

Definition. Let F = (T, $, 99, h) be a Friedrichs diagram and A C -Rr be 
a subset of the set Rt of lines of T and r be a map from Rt into 1R+. 

We define the quotient diagram Vat '■= {Ta, ^a, ^Af, hA) in the following 
way. To obtain the tree Ta we must tighten all lines from A into points. To 




(95) 



G(Or(r), ^?(/+(r)), ^7(/-(r)))(p)5(p - p') 




^P)'^9{{rJ-{r))) \P))- 



(96) 
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obtain we must remove all loops obtained by tightening all lines from A 
into the point. 

Now let us define ipAf- Joint all the vertices of T to A. We obtain a tree 
denoted by "^T. Let {C^T} be a set of all connected components of ^T. Let 
vq be a vertex of $a corresponding to the connected component C^T of ^T. 
Put by definition: 



X e '^^(^T')'^ "^^'^ . (97) 

Let us point out the notations in the previous formula. Rm is a set of all 
lines of $^ such that f'^ij-) and f~{r) are the vertices of C^T. {RT)r denotes 
the set of all lines of Rt such that the increasing path coming from f^{r) 
into f~^{r) contains r^. The symbol Vr denotes the set of all vertices v such 
that /+(r) >v> f~{r). hA{vo,r) = Yl h{v,r) + J2 '^tt- 

veV^Aj: r-reA; rTe(RT)r 

Definition. Let F be a Fricdrichs diagram. Let jFp be a space of all 
functions of external momenta of the diagram F of the form: 

V'(-Pext-), (98) 

where tp{...pext---) is a test function of external momenta. 

The convolution of the amplitude Ar{T){...pext---) with the function / e 
J^r we denote by Ar{T)[f]. 



8 The Bogoliubov — Parasiuk renormaliza- 
tion prescriptions 

Let for each Friedrichs diagram F = (T, $,(/?) A-p{T){...pext---) be some am- 
plitude. Fix some diagram F and let T' be some right subtree of the tree T 
corresponding to F. Let Tt' be a restriction of the diagram F on T' in obvious 
sense. Define the amplitude Ar^, ★ Ur{---Pext---) by the following formula: 
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^ Hie "--^ } n '^^(•••^^'■^-) 

reR' vev 

xAr^,{...p...). (99) 

In this formula V is a set of all vertices v such that v is not a vertex of Vt', 
i?' is a set of all Unes r of $r such that /"""(r) is not a vertex of T'. (i?^)^ 
is a set of all lines of T such that is not a line of T' and there exists 
an increasing path on T coming from f~{r) into f~^{r) such that this path 
contains r^. is a set of all vertices v of T such that v is not a vertex of T' 
and f~^{r) > v > f~{r). 

Let /lr('r)(p) be some amplitude. Put by definition: 

1 1 " 1 
ir(5i, := Ar{-, -){p) JJ -, (100) 

where n is a number of lines of Tp. Below we will consider the amplitudes 
Ar(s)[/] as distributions on (R+)" i.e. as an element of the space of tempered 
distributions S"((M+)"). Let ^{s} be a test function from S'((M+)"). The 
convolution of the amplitude Ar{s)[f] and the function ip{s) we denote by: 

:= J dsAr{s)[m^- (101) 

(R+)" 

The Bogoliubov — Parasiuk prescriptions. It will be clear below, 
that we can take into account only the diagrams F such that for each line 
tt of the corresponding tree of correlations T fji^rT ^ 3. Here Rrj. is a set 
of all lines r of F such that the increasing path on T which connects /~(r) 
and /"^(r) contains rx- Below we will consider only such diagrams. Other 
diagram can be simply subtracted by some counterterms At- 

According to the Bogoliubov — Parasiuk prescriptions we must to each 
diagram F (corresponding to the connected tree) assign the counterterm am- 
plitude Cr{s)[f] f G J-'r satisfying the following properties. 

a) (Locality.) C*r («)[/] is a finite linear combination of S functions cen- 
tered at zero and their derivatives. 

b) Let F be a Friedrichs diagram and T be a corresponding tree of corre- 
lations. Let A C Rrp and T' is some right subtree of T such that: 
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1) all lines of T such that is not a line of Rt' belong to A, 

2) All the root lines of T' do not belongs to A. 
Then 

Cr^Mif] = iCr'^„ ^ Ur)m], (102) 

where A' := ^4 n {Rt') and F' is a restriction of F on T'. 
c) 

0CAci?Tp 

where t — (ti,...,t„) = (■^,..., j-), the symbol C means here the strong 
inclusion and T is some subtract operator. 

d) The amplitudes Cr{s)[f] satisfy to the property of time-translation 
invariance, i.e. 

i Or{r)p'^t 

e ^ Crin, ...,Tn)[f] = Crin + t, ...,r„)[/]. (104) 

e) Let F be a Priedrichs diagram. Let 

Rrm]--=Urm]+ E CT^.ls)!/], and 
Rrm]--=Urm]+ E C'r,,(s)[/] + C'r(s)[/]. (105) 

The amplitudes Rr{s) are well defined distributions on (R"*")". 

f) The amplitudes Rr{s) satisfy the weak cluster property. This property 
means the following. Let f{.-.Pext---) be a test function. Then 

as a — > oo. Here pi is a projection of Pr to the x-axis. 
Put by definition for each diagram F: 

Ar(f) = E Ct,,, (107) 
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where ' in the sum means that all the root hnes of do not belong to A. 
Put 



At = ^ j dfAr{f){...pexf-)---at{Pext)---), 



(108) 



where the symbol F ~ T means that the sum is taken over all diagrams 
corresponding to T with suitable combinatoric factors. Suppose that the 
properties a) — f) are satisfied. Then At are the counterterms needed in the 
section . 

Theorem — Construction. It is possible to find such a subtract oper- 
ator T such that there exist counterterms C*r satisfying the properties a) — 



Note that it is not necessarily for us use not real counterterms. Indeed 
the evolution operator is real, so after renormalization we can simply take 
Re {RU){t, -oo)). 

9 Proof of the theorem-construction 

In this section we prove the theorem-construction from the previous section. 

Note that to prove our theorem we will use some ideas of the papers [10-12]. 
Before we prove our theorem let us prove the following 
Lemma 1. Let Li = S{R''), = <5((M+)"), k,n = 1,2,.... Let A{p) 

be some nonzero quadratic form on M^. Let T^^, i > be an one-parameter 

semigroup acting in Li drfined as follows: 



Let t > be some infinitely differentiable semigroup of continuous 
operators in L2. 

Let M be a subspace of finite codimension in L2. Suppose that M is 
invariant under the action of T^^, i.e. Vt > TfM C M. 

Suppose that there exist the linear independent vectors /i, .../; in L2 such 
that 



T/ : /(...p...) ^ e 



(109) 



Lin{{/„...,/;},M} = L„ 
MnLin{A,...,/J = o. 



(110) 
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and for each i = 1, t > 

T^h = h + ci,_J,_, + ... + ai/i + /, (111) 

for some coefficients aj_i, ai and the element / G M. 

Let g he a. functional on Li ® M such that g is continuous with respect to 
the topology on 5(M'=) x Suppose that V/ G Li ® M and Vt > 

{g, TJ) = {g, f) where T, = ® T^. 

Then, there exists an continuous extension g of g on S(R'') x 5'((]R'*')") 
such that V/ G Li (g) L2 and t> {g, Ttf) = {g, /). 

By definition we say that the functional h on Li® L2 is invariant if Vt > 
and V/gLi®L2 {h,TJ) = {hj). 

Proof of the lemma 1. At first we extend our functional g to the 
invariant functional ^ on Li ® L2 and then we prove that g is continuous. 

Let be a subspace of Li of all functions of the form, A{p)f{p), where 
f{p) is a test function. Let Mi = Lin{M U {/J}. Let 



/i:=(-T,2)Uo/i. (112) 
Let be a continuous functional on defined as follows: 

{k,ipip)Aip)) = -{g,ipip)(g)h). (113) 

Let k be an arbitrary continuous extension of k on whole space Li. The 
existence of such continuation follows from Malgrange's preparation theorem 
[13]. Now we define the continuous functional gi on Li ® Mi as follows: 

(^i,/®/i) = (fc,/), V/GLi. (114) 

According to (11131) we find that gi is an invariant extension of g on 
Li Ml. Step by step we can extend by the same procedure the func- 
tional g to the functionals ^2, ■■■■91 on Li ® Mi,...,Li ® M/ respectively, 
where M2 = Lin{MU {/„/,}},.. .,Mi = Lin{MU {/„/„..., /;}} respectively. 
Just constructed functional is separately continuous so it is continuous. The 
lemma is proved. 

Sketch of the proof of the theorem. We will prove the theorem by 
induction on the number of lines of the tree of correlations Tr corresponding 
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to the diagram T. It is evident that it is enough to consider only the diagrams 
with connected tree of correlations. 

The base of induction is evident. Suppose that the theorem is proved 
for all diagrams of order < n. (Order is a number of lines of the tree of 
correlations.) 

Let us give some definitions. Let ^{t) be a smooth function on [0, +cxd) 
such that < C,{t) < 1, C,{t) = 1 in some small neigborhood of zero and ^(t) = 

if t > Let us define a decomposition of unit {77^1(5) 1^4 C {1, ...,n}} by 
the formula 

Let ip{x) be some test function on real line such that ijjit) > 0, / ip{t)dt = 

1 and ■ip{t) = ii \t\ > jq. Put by definition: 

6x{x-X) = ^^^{^). (116) 

We have 

+00 

dX6x{x - X) = 1. (117) 



Let ^7v((M+)"), N = 1,2,..., be a subspace of 5((M+)") of all functions / 
such that / has a zero of order > at zero. Let "^{s) be a function of 
5((M+)"). 
We have: 

+00 

= J2 f / dsRr,,,{s\{i,...n}\A) 

Ac{l,...,n} { 



6,{l-\s\)^{Xs)r]A{^. (118) 

The inner integral in flllSp converges according to the inductive assump- 
tion. Therefore if ^I/(s) G S'7v((IR''')") and is large enough the integral at 
the right hand side of (11181) converges. So 
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defines a separately continuous functional on S'(R^-'^) ® S'Ar((M+)"'). / = / — !, 
where / is a number of external lines of F. To define a subtract operator T we 
must extend the functional (^r(s) [/], ^(s)) to the space S(R^^) (g) ^((M)") 
such that extended functional will satisfy to time-translation invariant prop- 
erty. To obtain this extension we use the lemma. In our case Li — S'(M^'^), 
L2 = A{p) = - Or{r)pl. is an operator acting in the 

5((M+)") as follows. 

Ttf{si, Sn) = /(j^^' ^2, Sn) if < ^, 

T,V(si,...,Sn)=0, ifs, >^. (120) 

The basis {/i, ...//} from the lemma is {s5"^...s^"770(s)}, mi, ...m„ = 1,2,3..., 
rui + m2 + ■■■ + rUn < N lexicographically ordered. We can now apply our 
lemma directly. 

Now let us prove the weak cluster property. Let p G R'^. Denote by 
P^i P'^ the projections of p to the x, y, z-axis respectively. To prove the 
weak cluster property it is enough to prove the following statement: for each 
connected diagram F the function {F-p{s){...pext---),'^{s)) defined by 

5(5^±Pe.t)(i^r(5)(...Pe.t-),^(5)) = {RT{^{...Pe.t...).^m (121) 

is a distribution of variables ...p1^^...p^^^... (constrained by momentum conser- 
vation law) which depends on ...p\^^... (constrained by momentum conserva- 
tion law) by the continuously differentiable way. We will prove this statement 
by induction on the number of lines of the corresponding tree of correlations. 
The base of induction is evident. Suppose that the statement is proved for 
all the trees of correlations such that the number of their lines < n. Let F 
be a diagram such that the number of the lines of the corresponding tree of 
correlations is equal to n. It is evident that if ^'(s) has a zero of enough high 
order at zero then (Fr(s)(---Pea;t---)> ^(^)) belongs to the required class (its 
enough to use our construction with decomposition of unit). Therefore we 
need to solve by induction the system of equations of the form: 

(i^±(Pe.t)')(Fr(«1(...Pe.t...),*(«1) 

= (Fr(s)(...fet...)47^t'*(«1)- (122) 
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According to Malgrangc's preparation theorem [13] we can choose the so- 
lution {FY{s){...pext---),'^{s)) such that it belongs to the required class if 
{Fr{s){...pext---), '^T^^{s)) belongs to the required class. Therefore the state- 
ment is proved. So our theorem is proved. 

10 Derivation of non ergodic property from 
main result 

Let us prove (more accurately as in introduction) that our system (Bose gas 
with weak pair interaction in thermodynamical limit) is non-ergodic system. 
Let us recall definition of ergodicity [15]. 

Definition. Consider a quantum system described by Hamiltonian H. 
This system is said to be ergodic if the spectrum of H is simple. 
This definition is equivalent to the following 

Definition. A quantum system described by Hamiltonian H is said to 
be ergodic if each bounded operator commuting with H is a. function of H. 

The generalization of this definition to the case when where exists some 
additional commuting first integrals is obvious [15]. 

It is to difficult to define a Hilbert space and Hamiltonian (as a self- 
adjoint operator in Hilbert space of states) of the system, in thermodynamical 
limit. So we give some new definition of ergodicity for this case which can be 
considered as some variant of last definition. Let us introduce some useful 
notations. Let V be an algebra of all Wick monomials with kernels from the 
Schwartz space, i.e F is a linear space of all expressions of the form 

/n m 
w{pi, ...,Pn\qi, ■■■,qm)Yla'^{Pi)dpiYla{qi)dqi, 
i=l 3=1 

w e S{R^^''+"'^), (123) 

where the multiplication is defined by canonical commutative relations. Let 
V be an algebraically dual of V. We say that the functional p E V is a, 
stationary functional ii^v E V p{[H,v]) — 0. Here if is a Hamiltonian of 

our Bose gas. We say that the functional p E V is a. translation-invariant 
functional if\fv E V p{[P,v]) = 0, where P is an operator of momentum of 
our system. Note that \/v E V [H,v], [P,v] E V. Denote by Vg the linear 
space of all Wick monomials of the form a[H,v] + 9 ■ [P,v]. Denote by 
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the space of all translation-invariant stationary states. We have V/ e V'{f e 
-v^ Vf G K f{v) = 0). Now let us introduce a notion of Gibbsian states. 
Let /3gM, /3>0,'yeM^. We define Gibbsian state on V formally by 

the following formula: 

^tr(ae-^(^+^^^))^ (124) 

where a eV, and Z is so called statistical sum: 

tr(e-^(W)). (125) 

This states corresponds to canonical distributions. Note that one of the 
basis statement of statistical mechanics states that there no difference which 
distribution we use: canonical or micro-canonical distribution. Bellow we 
will omit V, P at all formulas to simplify or notations. 

Let Vq be a subspace spanned by all Gibbsian states, i.e. is a set of 
all functionals {■) on V oi the form: 

(•) = E^"(-)/^- (126) 

a 

where the sum is understood in some generalized sense, for example it may 
be continuous (integral). It is evident that Vq (^Vs- 

Now we can give the definition of ergocity for Bose gas in thermodynam- 
ical limit. 

Definition. We say that our system is ergodic if each translation invari- 
ant stationary state can be represented as a superposition of Gibbsian states, 
i.e V: = V^. 

After these previous discussion let us start to prove our statement. Recall 
that we find non Gibbsian real stationary translation invariant functional (•) 
constructed as a formal power series on coupling constant A satisfying to the 
weak cluster property. This functional can be represented as follows: 

(•) = (-)o + (-)i + -, (127) 

where is a functional of zero order of coupling constant A, {■)o is a func- 
tional of first order of coupling constant A ant e.c.t.. 
Suppose that our system is ergodic. 

We do not suppose that the series {•)o + (•)! -|- ... converges, but we will 
work with it formally as with convergent series and find explicit formulas for 



33 



(•) under the assumption of ergodicity. Let us illustrate formal manipulation 
that we will use by several examples. 

oo 

Example 1. Let us calculate the sum '^x\ We do not suppose that 

i=0 

\x\ < 1. Denote by S the sum of this series. We have 

oo oo oo 

5 = ^a;^ = l + ^a;^ = l + x^a;* = l + x5. (128) 

1=0 i=l i=0 

Therefore 

S^l + xS, 
{l-x)S^l, 

s = (129) 
1 — X 

oo 

Example 2. Let us calculate the sum ^. Consider the function = 

i=l * ^ 

qqx + aix + a2x'^ + ... as a polynomial of infinite degree. Let us use the Viete 
theorem for this "polynomial". The roots of are Xi — in, i e Z \ {0}. 
According to the Viete theorem we have 

n ^i = C: 

iez\{o} 

E n ^^=0' 

j€Z\{o} iez\{o,j} 



E n 



i<j, i,jeZ\{0} k&\{0,i,j} 



(130) 



We find from these equations that 

E - = o. 



iez\{o} 



and 



i<j,ij6Z\{o} ' iez\{o} * jez\{o} ' 
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1/2 E 



1 



1 



6 



iez\{o} 



(131) 



Therefore 



6 



1 



(132) 



But Xn = Tm, so we finally have 



oo ^ 
i=l 



6 



(133) 



Such formal manipulation was widely used by Euler and others. Suppose 
that the set of such formal rules is enough large from one hand and does not 

contain a contradiction from other hand. These rules we call the Euler rules. 
If we can find the "sum" of some series by using the Euler rules then this 
series is called convergent in Euler sense. The "sum" of this series is called 
a sum in Euler sense. 

Let us prove that our functional (•) can be represented as follows (under 
the assumption of ergodicity): 



where are the "sums" of probably divergent series. The convergence (in 
the Euler sense) of this series will be proven below (under the assumption of 
ergodicy). The sum can be continuous (integral). 

Let {ea,a G 21} be a Hamele basis of Vs, Vs = Lin{eQ,,a; G 21}. Let 
P G ^}, 21 n 03 = be a completion of {ca, a G 21} to the Hamele basis 
of V, i.e. {ea, a G 21} U {ep, /3 G <B} be a Hamele basis ofV. V7 G 21 U 05 
let be an element of V such that f-yiey) — 1 and /y(ey) = if 7 7^ 7', 
7' G 2lU 03. 

An arbitrary functional p from V' now can be represented as a sum 



where are arbitrary numbers. Note that for arbitrary l^, l/s the right 

hand side of last equation is well defined because \/v & V fy{v) ^ only for 




(134) 



a 




(135) 
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finite number of elements 7 G 21 U *B. It is obvious now that an arbitrary 
element f & QVq (ergodicity) can be represented as follows: 

/ = EW/3, (136) 

where are arbitrary numbers. 

Vi = 0, 1, 2, ... we have the following representations 

But Vq; e 21 we have 

00 00 

i=0 1=0 

because Cq G Vs and (■) is a translation invariant stationary functional. 
Therefore Va G F 

(«) = E(-). = E(E 'afaia) + E 4//^(«)) 

i=o 1=0 aea /3e<B 

00 00 

= EE^«/"(«) + EE4//3(«) 

i=o aea i=o /3e® 

00 00 

-E(E^a)/"(«) + EE4/M«) 

00 

-E(E4//5)(«) (139) 

i=0 i3e^ 

in Euler sense. Finally 

00 

(«) = E(«)^' (140) 

where we put {■)[ — ^ s^/^ G Vq and our statement is proved. 

Let 1 be some enough large but finite subsystem of our system. Let 2 
be a subsystem obtained from 1 by translation on the vector / of sufficiently 
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large length parallel to the Let 12 be a union of the subsystems 1 

and 2. Let Ui, U2 and U12 be density matrices for the subsystems 1, 2 and 
12 respectively (which correspond to (■)). 

Let {ifil^, n = 1,2, ...} be a basis of eigenvectors of Hamilton operator for 
subsystem 1. Let{(p'^,m — 1,2,..} be a basis of eigenvectors of Hamilton 
operator for subsystem 2. Then the basis of eigenvectors of Hamiltonian for 
subsystem 12 is = V^n^V^m, Tn,n = 1,2,...}. According to the assump- 

tion that the systems 1,2 are enough large we find the following expression 

for density matrix of subsystems 1,2 in obvious notations. We have also the 
following expression for the density matrix of subsystem 12. 

Ca^®^- (142) 

But if / = +00 the weak cluster property implies that U\2 = Ui® 1/2- This 
leads to the following relations: 

a a p 

Here {En, n = 1,2,...} be a set of eigenvalues of Hi 2- But the set of 
sequences 

{e-fe} (144) 
is linear independent. So for each a we have 

e e 

According to the hnear independence of {e'"^} we find that Va 

Ca = 4. (146) 

So the series representing Cq, are convergent in the Euler sense and Va Cq, = 
0,1. 
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But we have — trC/ = 1. So for some (3 cp — 1 and Ca — H ck ^ (3. 

a 

We see that (•) is a canonical Gibbsian distribution 

(•) = (147) 

for some inverse temperature j3. 

Let us calculate now {a(k)a^ (k')) . We have constructed (•) by some 
Gauss state po described by some test function n{k). It follows from our 
construction of (•) that: 

{a{k)a+{k')) = n{k)5{k - k'). (148) 

But if momentum k is sufficiently large we can neglect by potential energy 
and find 

{a{k)a+{k')) = conste-^^5{k - k'). (149) 

If we chose n{k) such that n{k,) tends to zero as A; — > oo slowly than each 
Gauss function we obtain a contradiction. This contradiction proves non 
ergodic property of our system. 

Now let us discuss so called the Boltzmann ergodic hypothesis (1871). Let 
(■) be a translation invariant stationary functional on V such that Vt G M 
the functionals (e**^(-)e~**^) are well defined. The Boltzmann hypothesis 
states that for each such functional there exists an element (•)' G Vq such 
that Va e y 

^hrn^^|(e^*^ae-^*^)di=(a)'. (150) 



We see that according to Vg ^ Vq the Boltzmann hypothesis does not hold. 



11 Examples, chain diagrams 

In this section we consider by direct calculation some class of divergent di- 
agrams in Keldysh diagram technique. At first let us introduce the basis 
notion of the Keldysh diagram technique. 

Let us introduce the Green functions for the system 

p(r(*±(ii,a;i),...,*±(t„,x„))). (151) 
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Symbol H near means here that are Heizenberg operators. 

We require in nonequiUbrium diagram technique the following represen- 
tation for the Green functions 

Po(5~'T(*±(ti, xi), Xn)S)). (152) 

The symbol near means here that ^l'^ are operators in the Dirac rep- 
resentation (representation of interaction). The ^'-matrix has the form 

+00 

5-rexp(-i j V{t)dt), (153) 



and 



+00 

fexp(i J V{t)dt). (154) 



T is a symbol of the antichronological ordering here, po is some Gauss state 
defined by density function n{k) as usual. 

Let us recall the basic elements of nonequilibrium diagram technique. 
The vertices coming from T-exponent are marked by symbol — . The vertices 
coming from T-exponent are marked by symbol -|-. There exist four types of 
propagators 

Gfiti - t2, xi - X2) = po(^(^i, a;i)^+(t2, X2)), 

Go^ih - t2, Xi - X2) = Po(^^(t2, X2)^(tl, Xi)), 

Go"(^i -^2,a;i -0:2) =po(T(^(ti,a;i)^+(t2,a:2))), 
G^+ih - t2, xi - X2) = Po(T(^(ti, a;i)^+(t2, X2))). 



(155) 



Let us write the table of propagators 

^o^^t^^) = J ^(27r)5(a;-a;(A:))n(A:)e-^("), 
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{27iy^u -uj{k) +iO uj-Lo{k)-iO 
d^k n{k) l + n{k) 



{27r)^^LU -Lu{k) + iO lu - Lu{k) - iO 



d^k l + n{k) n{k) 



(156) 



■i{u>t—kx) 



11.1 D ivergences 



fig. 1 



Or 



• • • 



A typical example of divergent diagram is pictured at figl. 
The ovals represent the sum of one-particle irreducible diagrams. These di- 
agrams are called chain diagrams. Let us suppose that all divergences of 
self-energy parts (ovals) are subtracted. The divergences arise from the fact 
that singular supports of propagators coincide. At first we consider diagrams 
with one self-energy insertion (one-chain diagram). These diagrams are pic- 
tured at fig. 2. 



These diagrams are analogous to 

one-loop diagrams in quantum field theory. 

The aim of this section is to prove that the Green functions can be made 
finite by the following renormalization of the asymptotical state: 



fig. 2 




Po(-) ^ ^Po(e 



/ h(t)dt 



(•)), 



(157) 



— oo 



where 



h 
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h{k) is a real- valued function and 

+ 00 

- / h{t)dt 

Z^poie - ). (158) 

11.2 Proof of the existence of divergences in the theory 

Suppose that there are no divergences in the Keldysh diagram technique if 
''T'i^) — for positive a, /3. Therefore the Green function 

p(S-\S^t(ti, xi)*o(i2, X2))) (159) 
is translation invariant. So the density matrix 

Pt{xi,X2) := p{S-\S^+{t, xi)^o{t, X2))) (160) 
is an integral of motion. Let 

pt{k)^ J d'xpt{0,x)e"'^. (161) 

In zero order of perturbation theory p[k) = n{k). But if there are no diver- 
gences in Keldysh diagram technique it is possible (see [14]) to derive the 
following kinetic equation for p{k) 



I 



dptjk) 
dt 

«^(P,Pl|P2,P3){(l + P(P))(1 + P{P1))P{P2)P{P3) 

-p{p)p{pi){l + P(P2))(1 + P(P3))}. (162) 
The right hand side of this equation is equal to zero only if 

Pik) = (163) 

e«^+^ + 1 

for some a, P (a > 0, /3 > 0). But n{k) = p{k) in zero order of perturbation 
theory, so n{k) has a Bose- Einstein form. This contradiction proves our 
statement. 
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11.3 Regularization 

Let us now introduce regularization. Note that 



1 



X 



— TT 



X + ie x'^ + e 
Therefore we use the following regularization 

S{uj — uj{k)) — > — 
1 



(164) 



n 



--: 5s{uj - uj{k)), 
1 



G^^{t2 — ti,X2 — Xi), 
■ G^'{t2-ti,X2-Xi). 



'uj — uj{k) {uj — uj{k)y + e"^) ' oo — uj{k) 

11 A Some simple relation on the Green functions 

Lemma 1. The following equalities hold 

G {ti - t2,Xi - X2) 
G^ [ti — t2, Xi — X2] 

Proof. We have 

G—{ti - t2, Xi - X2)* = Po(^(^//(tl, Xi)^Hit2, X2))) 
= po(T(^+ (ti, Xi)^H(t2, X2))) = G'++(t2 - tl, X2 - Xi) 

So the equality 11661 is proved. We have 

G+- {ti - t2,Xi - X2y = Pq{:^ H{ti,Xi)^+(i2,X2)y 
= Poi'^Hit2,X2)^Hitl,Xi)) = G+-{t2-ti,X2 - Xi). 

So the equality 11671 is proved. 
The Lemma is proved. 
It is easy to prove the following 
Lemma 2. The following equality holds 

x) = e{t)G—{t, x) + e{-t)G++{t, x), 
G-+{t, x) = e{t)G++{t, x) + e{-t)G—{t, x). 

Let us introduce the following matrix 



(165) 



(166) 
(167) 



(168) 



(169) 



(170) 



G 



G-' 



G- 



(171) 
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Let us introduce the similar matrix for the self-energy operator 



Dyson equations in Fourier representation have the form 

G = Gq -\- GqYjG. 

We have from these equations that 

S = Go ^ - C-\ 

or in the matrix form 



detGo 



G, 




Gn 



-Gq 

Gn 



1 



detG 



G- 
-G- 



G+ 



It follows from Lemma 1 that 



G++{Lo,p)=G-iuj,py 
G+-{uj,p) = G+-{uj,p)\ 
G-+{co,p) = G-+iuj,py. 

Therefore detGo, detG are real and we have the following lemma. 
Lemma 3. 



S++(t2 - ti,X2 - Xi), 
- S+"(t2 - ti,X2 - Xi). 



s (ti - t2,xi - X2y 
- 12, xi - X2y 

The following Lemma holds. 
Lemma 4. 



S++(cu,p) + E-{uj,p) = -S-+(cu,p) - E+-(cu,p). 

Proof. The statement of lemma follows from the Dyson equat 
and the following two obvious equalities: 

G++(a;,p) + G-iLu,p) = G-+(cu,j9) + G+-(cu,p), 
G++(a;,p) + Go-(^,p) = Go+(^,p) + G+-{uj,p). 
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11.5 Calculation of the propagators in one-chain ap- 
proximation 

Lemma 5. The following limit equality holds (in the sense of distributions): 

\im{-Sfr{x) = reg, 

lim{- ^ - -S{x)} = reg, 

\mi5e{x)Ve{-) = reg, 

xeR. (181) 

Here reg means some correct distribution. 

Proof. Let f{x) be some test function with compact support. We have 



= ^ / + ^/'(O) + ^''{'i^)}''^ (182) 

for some smooth bounded function ip{x). We have 

/ '^'W(^) = ^ / ^^^^{/(0) + 5a;/'(0)+£V^(.a;)} 

{-Jj^^dx}m + 0{e). (183) 



TT 

But 



2 ' 



dx = -. (184) 



So 



(X2 + 1)2 2 



5'(^)/(^) = ^/(0) + O(^). (185) 



So the first equality is proved. One can prove other three equalities in the 
same way. 
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Therefore wc sec from the Lemmas 1,2. that we can consider only the 
function G {t^x). But the function G {t^x) can be represented as a sum 
of chain diagrams. At first let us consider the diagrams with one self-energy 
insertion (one-chain diagram). We have G~~ = ^ H^^ , where the diagrams 

for Hl^ are presented at the fig. 2. We have the following representation for 
the divergent parts of these diagrams. 

= 2.E-(c.,pMp)(l + n(p))i%-a;(p)) 

+27rE++(w,p)n(p)(l + nip))-S{uj - uj{p)). (186) 

e 

We see that the divergent part of these two diagrams is real (because E = 

(E++)*). 

Let us consider the singular part of other two diagrams presented at fig 

3. 

fig. 3 




We have 

{H~^)div{uj,p) + {H^~)div{uj,p) 
- 7r(27r)(l + 2n(p))(l + n{p))S^,{uj - a;(p))E-+(a;,p) 
+7r(27r)(l + 2n{p))n{p)5^^{uj - u{p))E+-{uj,p) 

= 7r(l + 2n(p))(l + n{p))-6{uj - uj{p))T.-+{uj,p) 

s 

+.(1 + 2nW)n(p)ii(. - .(p))E-(.,p) + 0(.). (187) 

We see that {H--)diy{uJ,p) + {H++)diy{uJ,p), {H-+)div{uJ:P) + {H+-)diy{uJ,p) 
are real. 

We will use the dotted line for lines which connect creation- annihilation 
operators with operators arising from the vertex: J h{k)a^ {k)a{k)(Pk (see 
fig- 4). 
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can be subtracted by the following counterterm: 

;i + 2n(p)) 



So the divergences in = ^ Hl^ 



%) = E++(a;,p) + E--(u;,p) + 



2n{p){l + n{p)) 
x{(l + n(p))S-+(a;,p) + n(p)S+-(cc;,p)}. 



By using Lemma 4 we have: 



{(l + n(p))E-+(u;,p) -n(p)E+-(a;,p)}. (189) 

The left hand side of this equation can be rewritten as follows (in approxi- 
mation used in [14]) 

h(p) = _1±^!^M— -5i(p), (190) 
2n(p)(l + n(p)) ^ ^ 

where St{p) is a scattering integral. So h{p) ^ for non-equilibrium matter. 

Analogously one can consider two-chain diagrams presented at fig. 5 by 
direct calculation and prove that the divergences can be subtracted by the 

fig. 5 



counterterms of the asymptotical state see cite. 

Note that there arises the phenomenum of overlapping divergences in this 
example. 
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12 Notes on Bogoliubov derivation of Boltz- 
mann equations 

In this section we study the problem of boundary conditions in Bogohubov 
derivation of kinetic equations [6]. Let us consider N particles in M'^. Let 
Qi be a coordinates of particle number i, and Pi be a momentum of particle 
number i, i — 1, N. Suppose that particles interact by means of the pair 
potential V{qi — qj). We suppose that V belongs to the Schwartz space. Let 
Xi = {Pi, Qi) be a point in the phase space F. Let f{xi, a;„) be a distribution 
function of N particles. If wc want to point out that f{xi, ...,Xn) depends 
on t we will write f{xi, ...,XN\t). Let 




(191) 



and 



/2(xi,X2)= J dX3, ...,dxNf{xi, ...,Xn) 

be marginal distribution functions. Put by definition 

pi{xi) = Nfi{xi), 

P2{Xi,X2) = N'^f2{Xi,X2). 

li A is a, function on the phase space F, F = M.^^ and 



(192) 



(193) 



N 




(194) 



i=l 



then 




(195) 



Now ii A is a function on the phase space 




(196) 
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in the limit of large N we find 

{A) ^ J dxidx2P2{xi,X2)A{xi,X2). (197) 
Let us introduce also three-particle distribution function: 

Mxi,X2,X3) = / dX4...dXNf{xi,...,XN)- (198) 



Let us derive the equation for f{x). At first let us write equation of motion 
for f{xi, ...,xn). We have 

o n 

— /(xi, XN\t) + } — Vi/(a;i, Xn\t) 

at ^-^ m 

1=1 

dV{qi - qj) dfjxi, ...,Xn\t) 

This equation is only an infinitesimal form of the Liouville theorem. Let us 
multiply this equation by and integrate over dx2, ...,dxN- Suppose that 
/(xi, ...,xn) is a function of rapid decay of momenta. This assumption allow 
integrate over pi by parts. We find: 

d p 

— pi(xi|i) + —Vpi{xi,t) 

ot m 



+ I dX2 — ^P2{Xl,X2\t) 

J mdq2 



^^m<ll-<l2)dp2{x,,X2\t)_ (200) 

oqi dpi 

Note that we have kept here boundary term. Let us now talk about deriva- 
tion of kinetic equation. According to the standard prescription we put 
P3{xi,X2,X3) = in equation for p2{xi,X2)- We find the following equation 
for P2: 

j^P2{xi{t),X2{m = Q, (201) 
where {xi{t),X2{t)) is a solution of corresponding two-body problem. 
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Condition of correlation breaking. We consider only translation- 
invariant matter in purpose of simplicity. Usual correlation-breaking condi- 
tion has the form 

P2(a;i,X2|0) = h{p[{xi,X2))h{p2{xi,X2)). (202) 

Here h is a function on momenta-space of one particle. We consider only 
translation-invariant gas, so h depends only of momentum. 

p'i{xi,X2) and p'2{xi,X2) are momenta of particles 1 and 2 at t = — oo if 
at if: = their coordinates and momenta was Xi and X2 respectively. 

Proposition. 

^P2(xi,X2) = 0. (203) 

Indeed, according to (12011) 

P2(Xl,X2|t) =P2(X?,X°|0), (204) 

where and X2 are phase coordinates of particles 1 and 2 respectively at a 
moment t = 0. Therefore 

P2(xi,X2|t) = /i(Pi(x?,x°))/i(p'2(x?,x°)). (205) 

But the points and a;^ come to the points xi and X2 after the time t. So 

{p'i{xi,x°2),P2{x1,x9,)) = {p[{xi,X2),P2{xi,X2)), and 

p2{xi,X2\t) = h{p[{x^^, xl))h{p2{x[, xD) = 
h{p[{xi,X2))h{p2{xi,X2)) = P2(xi,a;2|0). (206) 

In result 

P2(a;i,X2|t) = p2(a;i,X2|0). (207) 

The proposition is proved. 

It follows from equations (I20ip and (12031) that: 

(-Vi + ^V2)/2(a;i,X2|t) 
m m 

dV{qi-q2) d dV{qi - ^2) d 

= + -7^)f2{Xi,X2\t). 208 

dqi dpi dq2 dp2 

49 



The function h{p) can be found from the following equation 

pi{x) = lim ^ / p2{xi,X2)dx2. (209) 

N-^oo iV J 

But in zero order of gas parameter the particles are free and 

pi(x) = h{x). (210) 

Formula (12081) is usually used for transformation of r.h.s. of equation (12001) 
to the scattering integral. From other hand the equation 

^P2(Xl,X2) =0 (211) 

shows that there is no irreversible evolution in the system. From other point 
of view we will show that the last term in the left hand side of (12001) is equal 
to the scattering integral. 

For simplicity we will show the case vi = 0, v = ^. The general case 
can be reduced to this case by means of Galilei transformation. So let us 
consider the integral 

/= lim [ d^p2 [ rfg2-^P2(0,0,p2,g2), (212) 
R^oo J J mdq2 

Vr 

where Vr is a ball of radius R with the center at zero. Let us integrate over 
dq2 by using Gauss theorem. We find 

/= lim / d^p2 f dS^ cos i:p2i0,0,p2,q2). (213) 

Sr 

Here Sr is a boundary of Vr and is an angle between two rays: first of 
them is parallel to p2, second starts from zero and passes throw q2. We have 

d p2 dS — cos ip X 
J 

Sr 

h{p[{0,0))h{p'2{p2,q2)). (214) 

Let us suppose that the particles scatter only then they are not too far from 
each other. Then 

h{p[{0,0))h{p'2{p2,q2)) = Pi(p2)pi(0) (215) 
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for all q2 G Sr\0, where C is a small neighborhood of the point go '■= ^ 
Sr. Diameter of O is approximately equal to diameter of suppF. Therefore 
the integral / is not equal to zero and equal to 



where b := q2 — %■ But the right hand side of (12161) is a usual scattering 
integral. 

Therefore if we keep boundary terms in BBGKI-chain we obtain the ki- 
netic equations without scattering integral. 

13 Conclusion 

In the present paper we have developed the general theory of the renormal- 
ization of nonequlibrium diagram technique. To study this problem we have 
used some ideas of the theory of /^-operation developed by N.N. Bogoliubov 
and O.S. Parasiuk. 

We illustrate our ideas by simple example of one- and two-chain diagrams 
in Keldysh diagram technique. 

We want to illustrate in this paper the following general thesis: to prove 
that the system tends to the thermal equilibrium one should take into ac- 
count its behavior on its boundary. In the last section we have shown that 
some boundary terms in BBGKI-chain which are usually neglected in Bogoli- 
ubov derivation of kinetic equation compensate scattering integral in kinetic 
equation. 
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